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Abstract 

Using the relation established by Johnson-Zumbrun between Hill's method of aprox- 
imating spectra of periodic-coefficient ordinary differential operators and a generalized 
periodic Evans function given by the 2-modified characteristic Fredholm determinant of 
an associated Birman-Schwinger system, together with a Volterra integral computation 
introduced by Gesztesy-Makarov, we give an explicit connection between the generalized 
Birman-Schwinger-type periodic Evans function and the standard Jost function-type 
periodic Evans function defined by Gardner in terms of the fundamental solution of 
the eigenvalue equation written as a first-order system. Fhis extends to a large family 
of operators the results of Gesztesy-Makarov for scalar Schrodinger operators and of 
Gardner for vector-valued second-order elliptic operators, in particular recovering by 
independent argument the fundamental result of Gardner that the zeros of the Evans 
function agree in location and (algebraic) multiplicity with the periodic eigenvalues of 
the associated operator. 

1 Introduction 

The purpose of this note, generalizing results of [GM04, GLMZ05, GLM07J for asymp- 
totically constant-coefficient operators and of Gesztesy-Makarov [GM04] for of periodic- 
coefficient Schrodinger operators, is to give an explicit connection between two objects 
related to the spectra of periodic-coefficient operators on the line, namely, a Birman- 
Schwinger-type characteristic Fredholm determinant -D(A) introduced in \JZ\ |BJZJ and the 
Jost function-type periodic Evans function E(X) of Gardner [G] . 

Both of these objects are analytic functions whose zeros have been shown in various 
contexts to agree in location and multiplicity with the eigenvalues of the associated operator. 
Thus, on the mutual domain for which these properties have been established, it follows 
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that they must agree up to a nonvanishing analytic factor. Here, we determine explicitly 
this nonvanishing factor, thus illuminating the relation between the two functions while 
extending the key property of agreement with eigenvalues to the union of the domains 
on which it has been established for each function separately. Knowledge of this scaling 
factor is useful also for numerical approximation of spectra as in [BJZ|. IBJNRZH IB JNRZ2 . 
allowing convenient comparison of different methods. 

Our method of proof proceeds by comparison with an interpolating first-order Birman- 
Schwinger determinant, using a simple version of the Volterra integral computation of 
[GM04] together with the relation established in [JZJ between Birman-Schwinger determi- 
nants and Hill's method, a spectral Galerkin algorithm for finite approximation of spectra. 
We carry out the analysis here for the simplest interesting case of a second-order elliptic 
operator; however, our arguments extend readily to the full class of operators discussed in 

Consider the eigenvalue problem for a general second-order periodic-coefficient ordinary 
differential operator L = Bq^x) -1 [d 2 + d x A\{x) + Aq(x)), written in form 

(1.1) {d 2 x + d x A x {x) + A (x) - XB (x))U = 0, 

where L is defined on complex vector- valued functions U G L 2 [0, X] with periodic boundary 
conditions, Aj,Bj G L 2 are matrix- valued and periodic on x E [0,X], and B$ is positive or 
negative definite in the sense that its symmetric part 'RBq := (1/2) (Bo + Bq) is positive or 
negative definite0 (Note that, by Lyapunov's Lemma, definiteness can always be achieved 
by a variable-coefficient change of coordinates, provided that the eigenvalues of Bq are of 
purely positive or purely negative real part.) Taking the Fourier transform, we may express 
(jl.ip . equivalently, as an infinite-dimensional matrix system 

(1.2) {V 2 + VA + B)U = 0, 

Bjk = - A-Boj-fc, Ajk = M,j-k, = (2ir/X)5jij, Uj = U(j), where / denotes the 

discrete Fourier transform of /. (Here and elsewhere i = \f— 1.) We shall alternate between 
these two representations as is convenient. 

2 The Birman— Schwinger Evans function 

By a Birman-Schwinger-type procedure, applying (d 2 — 1) _1 to (jl.ip on the left, we obtain 
an equivalent problem 

(2.1) (J + K{\))U = 0, 

where K = K l + K , with K\ = d x (d 2 - l)" 1 ^, K = (d 2 - 1) _1 (A) + 1 - XB ). 
Lemma 2.1 (|JZj). For Aj,Bj G L 2 , the operator K is Hilbert- Schmidt. 

The assumption of symmetry of Bo made in [JZ] is not necessary for the arguments there. 
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Proof. f [JZ| ) Expressing K in matrix form /C with respect to the infinite-dimensional Fourier 
basis {e 2wl ^ x ^ x /y/X}, we find that K,\jk = 1 +{2n/x)^j 2 ^ 1 > J ~ fc 1 anc ^ srmuarr y f° r where 
denotes Fourier transform and i := v— 1. Taking without loss of generality X = 27r, we 
find that 

J 2 J 2 

(2.2) Uplift = H/Cillsa =Z] (l + j2)2 Zj Al (j~ fc )| 2 = (l +J -2) 2 H A llU 2 W < +°°' 

and similarly for /C 2 . □ 
Definition 2.2. We define the Birman-Schwinger Evans function as 

(2.3) D(X) :=det 2 (I-if(A)), 

where det 2 denotes the 2-modified Fredholm determinant, defined for Hilbert-Schmidt per- 
turbations of the identity; see Appendix [A] for a review of the relevant theory. 

Proposition 2.3 iBJZj). For Aj,Bj G L 2 , D is analytic in X; moreover, its zeros 
correspond in location and multiplicity with the eigenvalues of L := B^ 1 (d 2 + d x A\ + Aq). 

Proof. Analyticity follows by analyticity of the finite-dimensional Galerkin approximations 
by which the Fredholm determinant is defined, plus uniform convergence of the Galerkin 
approximations, a consequence of (IA.3j) , Correspondence in location is immediate from the 
fact that det 2 (/ — K) = if and only if (I — K) has a kernel, while correspondence in 
multiplicity may be deduced by consideration of a special sequence of Galerkin approxima- 
tions on successive eigenspaces of L, for which the Galerkin approximant can be seen to 
be a nonvanishing multiple of the characteristic polynomial for the restriction of L to these 
finite-dimensional invariant subspaces. See |JZj for further details. □ 



3 Hill's method 

Hill's method consists of truncating (|1.2j) at wave number J, i.e., considering the (2j + 1)- 
dimensional minor \j\ < J, and solving the resulting finite-dimensional system to obtain 
approximate eigenvalues for L: that is, the eigenvalues of the (2 J + 1) x (2J + 1) matrix 

(3.1) Lj = (13 ,j)- 1 (Vj + VjAj + Bj), 

F>o,J,jk = -Bo,j-fc; where subscripts J indicate truncation at wave number J, or, equivalently, 
restriction to the Jth centered minor as described above. 

Left-multiplying by (T>j — 1) , similarly as above, we may rewrite the truncated eigen- 
value equation (Lj — X)U = (V 2 + VjAj + Bj)U = equivalently as (I + K, j)U = 0, where 
JCj = K\ t j + K.2,j is the truncation of the Fourier representation K. = K,\ + IC2 of operator 
K, that is, Ki = Vj(V 2 - /j-Uy, K 2 = (V 2 - I)~ l {Ao,j + 1 - XB 0J ). We define the 
truncated Birman-Schwinger Evans function, accordingly, as 



(3.2) Dj(X) := det 2 (I-£j) = det(Vj - 1) det £ ,je tr/Cj(A) det(Lj - A). 
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Proposition 3.1 ( |JZj ). Each Dj is analytic in X; moreover, the zeros of Dj correspond 
in location and multiplicity with those of Lj. 

Proof. Immediate, from (|3.2p . the assumed uniform positive definiteness of Bq, hence also 
of Bo and £>o,J> and properties of the characteristic polynomial. □ 

Proposition 3.2 ([JZ]). ForAj,Bj G L 2 , Dj(X) -> D(X) as J -> oo, uniformly on \X\ < R. 

Proof. By diagonality of P, Dj.Aj = (2XA)j, hence Dj is a sequence of Galerkin approxima- 
tions of the Fredholm determinant D = det2(4 — K) on a complete set of successively larger 
subapaces. By definition of the Fredholm determinant, therefore, Dj — > D as J — > oo. □ 

Corollary 3.3 (Convergence of Hill's method |JZ| ). For Aj,Bj G L 2 , the eigenvalues of Lj 
approach the eigenvalues of L in location and multiplicity as J — > oo, uniformly on \X\ < R. 

Proof. Immediate, from Lemma 12.31 Propositions I3.l l and 13.21 and properties of uniformly 
convergent analytic functions. □ 

4 The Evans function of Gardner 

Rewriting (|l.ip first-order system 

(4.1) B,W.mw. := ((AB -A - d x A{) - A ) • 
we may alternatively define an Evans function by shooting, as done by Gardner. 
Definition 4.1. We define the Jost-type Evans function of Gardner, following [Gr) . as 

(4.2) E(X) :=det(*(Z) -/), 

where denotes the fundamental solution of (j4.ip . satisfying ^(0) = /, and X is the period. 

Proposition 4.2 ([2]). For Aj,Bj G C , D zs analytic in X; moreover, its zeros correspond 
in location and multiplicity with the (periodic) eigenvalues of L := B ~ 1 (d 2 + d x A\ + ^4o)- 

Proof. Analyticity and agreement in location follow immediately by analytic dependence 
of solutions of ODE and the fact that A is an eigenvalue precisely if there exists a solution 
W(x) = ^{x)Wq of (j4TT|) for which W(X) = ^(X)W = W . Agreement in multiplicity 
may be shown by a more detailed argument based on choice of a Jordan basis as in jjSj. □ 



5 CONNECTIONS 



5 



5 Connections 

Define now the intermediate first- order Birman-Schwinger Evans function 

(5.1) F(X) := det 2 ((d x - l)" 1 ^ - A)) = det 2 (I - K), 

where K := {d x - 1) _1 ( A - !)• 
Theorem 5.1. For Aj,Bj £ L 2 , 

(5.2) F{\) = 7 (e* - l)~ 2n E(X), 
where 

( 5 _3) . = e IZ^r J,f (tr(A)-2«)(y)dy = e ^(tT(A 1 , ave )+2n)X ^ 



and Ai^ ave denotes the mean over one period of A 



l- 



Remark 5.2. From ([5jZ])-([0]), we obtain F(A) ~ E(X)e~ £ tr ( A )(^ = det(J - ^{X)~ l ) 
as X — > oo, so that F is asymptotic to a "backward" version of Gardner's Evans function. 

Proof. By a straightforward computation taking into account periodicity together with the 
jump condition G(y + ,y) — G(y_,y) = 1 at x = y, we find that the Green kernel associated 
with (d x — 1) _1 is semiseparable, of form 



(5.4) G(x,y) 



fzpr, x > y, 

X + x-y 

-jzp^, x<y. 



Following [GM04], we may thus decompose G as the sum G(x,y) = f(x)g(y) + J(x,y) 
of a separable kernel 



(5-5) f(x)g(y): f(x) = e x , g(y) 



e x e~y 
l-e x 

and a Volterra kernel 



(5.6) J(x,y) 



e x y , x > y, 
0, x < y. 



Likewise, the kernel k(x, y) of IK decomposes as 

(5.7) y) = f(x)g(y)(K(y) - 1) + H(x, y), 

where H(x,y = J(x,y)(A(y) — 1) is again a Volterra kernel, vanishing for x < y. 
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Using the Fredholm determinant facts (|A.5|) . (1A,6|) compiled in Appendix [Aj together 
with the fact that det 2 (I — H) = 1 for a Volterra operator H$ we find (here, with slight 
abuse of notation, denoting operators by their associated kernels) that 

det 2 (l-(d x -l)-\A-l)) 

= dei 2 (l-H-fg(A-l)) 

(5.8) = det 2 (/ - F)det 2 (j - (J — H^fgiA - i)) e -tr(tf(J-*r l /*(A-i)) 

= det 2 (/- (I-H^fgiA- 1) ) e tr(/< ? (A-i)) e -tr((/-H)-V«/(A-i)) 

= 7 det c2 „ (/ - ( 5 (A - 1), (/ - H)- l f)) , 

where (•,-) denotes L 2 [0,X] inner product and 7 := e^ r (^ A ~ 1 )'-^) = e eX - 1 *■ t r (^i,<")e)+2n)x _ 
Computing 

(I - H)*(a;) = *(z) - / e :E ~ s/ (A- l)*(y)dy, 

Jo 

and recalling that (3a; — 1)\I/ = (A — l)^, where denotes the fundamental solution of 
(d x — A)^ = 0, ^(0) = /, so that, by Duhamel's principle/ variation of constants, ^(x) = 
e x + /* e x - y (A-l)^(y)dy, we find that (I-H)^ = e x = f, or (I-H)~ 1 f = Substituting 
into (|5.8p . we obtain, noting that := e~ y ^f(y) satisfies (9^6 = (A — 1)6, 



1 


-e x 




e x 






1 


-e x 




e x 






1 


-e x 







det 2 ( J — (d x - l)~ x (A - 1)) = 7det C 2« [I — { Y^TJc ) L {^-^K v ^(y)dy 

(5.9) 



7 det c2 „(l-( T ^ TF )e-^(y)|^ 



= 7 (e x -l)- 2n det c2 n(M/(X)-/), 
yielding the result. □ 

Now, define the truncated first-order Birman-Schwinger Evans function as 
(5.10) Fj(X) :=det 2 (/-/Cj), 

where Kj denotes the Jth Fourier truncation of the Fourier transform fC of K.. 

Lemma 5.3. For Aj,Bj G L 2 ; -Fj(A) — > F(X) as J — >■ oo ; uniformly on |A| < i?. 

Proof. As in the argument of Proposition 13.21 this follows essentially by the definition of 
the Fredholm determinant. □ 



Heuristically a strictly lower triangular perturbation of the identity; H has no nonzero eigenvalues. 
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Lemma 5.4. Similarly as in (|3.2p . we have 

(5.11) Fj{\) = det(Pj-l) 2 deti3 ,je tr,ej Wdet(Lj-A). 

Proo/. Using (J-K) = (9 X - 1)~ 1 (9 X - A), A(A) = _ ^ _ _^ ) to obtain 



Vj-Ij 1 



Vj -Ij \ 



and, using the block determinant formula det ^ = det(ca + 6), 

det ( Bj + V (k^h v~I j a) = det + ^ + ^ + B ^ 

= det (Vj + Pj^lj + Bj) 
= det B , j det (L j - A). 



Corollary 5.5. For Aj,Bj € L 2 ; F(X) = ee S ~ 5 D(X), w/iere ; /or X = 2vr, 

5:= tr(i 1 , + 2/)(l+^-^-), 

i<b'l 

(5.12) 6:=-tr(A lfi + I-\B ,o)Y^p-ri> 

j 

e := lim TT (l + 



Proof. We may readily calculate that 

tr/C = tr(ii j0 ) Yl _ l_ i + tr (^o + 7 " A ^o,o) _ =n i 

(5.13) 

= -tr(i li0 + /-APo,o) £ -2-y 



□ 



V + 
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and 



tr/C 



-(Vj-Ij)- 1 I 

(Pj-Ij)" 1 ^!-/) 



(5.14) 



tr(i lj0 + 21) V T — • 

\3\<J J 

tr(i 1)0 + 2/)(-l+ £ (tt^-t-- 



in — L in 

i<\j\<J 



tr(ix,o + 2/)(l+ 53 t4 



7 + ^7 

both uniformly convergent by absolute convergence of ^ 1/j 2 - Similarly, 

fa'" 1 ) 2 = TT fi 4- „ e E„<^+oO- 2 ) „ e E b -|<, o(i" 2 ) 



lil<J b'l<J 



may be seen to be uniformly convergent. Comparing (|3.2p and ()5.11|) and taking the limit 
as J — > oo, we obtain the result. □ 

Collecting information, we have our final result, relating D, E, and F. 

Corollary 5.6. For Aj,Bj G L 2 , and X = 2ir, 



(5.15) D(X) = —F(X) = —^e* - l)~ 2n E(X), 

e e 



where 7 is as in (|5.3p and 5, 5, and e are as in (|5.12p . 

Remark 5.7. Since ^—^^i(e x — l)~ 2n is analytic in A and nonvanishing, Corollary 15.6 
together with Proposition 12.31 gives an alternative proof of Proposition 14.21 



A 2-modified Fredholm determinants 

In this appendix, we recall for completeness the basic properties of 2-modified Fredholm 
determinants of Hilbert-Schmidt perturbations of the identity [CCK971 ICGKOOl IGK691 
ISi77l IST05] . For a Hilbert space H, the Hilbert-Schmidt class &2C%) i s defined as the set of 
linear operators on % for which 

WMW-H) ■= J3l(^,e fc >| 2 = tr(AM) < +00, 
where {ej} is any orthonormal basis. Evidently, || • ||e 2 (^) is independent of the basis. 
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For a finite-rank operator A, the 2-modified Fredholm determinant is defined as 
(A.l) det 2>n (I n - A) := det n ((I n - A)e A ) = det H (I H - A) e tT ^ A \ 

where det% denotes the usual determinant restricted to Range(A). There holds 
(A.2) e - c \\ A \\l 2 m < |det2, w (/« - A)\ < e cm2 ^m 

and 

(A.3) |det 2 ,«(/w - A) - det 2 ,w(/« - B)\ < \\A - £|| B2(w) e c[l|A|lB 2(w)+ll B llB 2 («) +1 ] 2 , 

where C > is a constant independent of the dimension of the space. 

For A G B 2 (T-L), the 2-modified Fredholm determinant is defined as the limit 

(A.4) det 2 , n (lH ~ A) := lim det 2 , Wj (/« 7 - Aj), 

J— >oo 

where Hj is any increasing sequence of finite-dimensional subspaces filling up T~L, and Aj 
denotes the Galerkin approximation Py^jAly^j, where Pj is the orthogonal projection onto 
Hj. Equivalently, det2,n(^'H ~ A) := 11^(1 — aj)e aj , where aj are the eigenvalues of A. For 
A € B2, (I - A) is invertible if and only if det 2 (I -A) / 0. 

Properties (|A.1[) — (|A.2[) are evidently inherited by continuity. Likewise, from the corre- 
sponding properties of finite-dimensional determinants, we obtain in the limit 

(A.5) det 2 , H ((I-H ~ A)(I H - B)) = det 2>n (I-H ~ A)dei 2>n (In - B) e - tr «( AB ) 

for all A, B G B2 (here we are implicitly using the fact that AB is in trace class, with 
||^4-B||bi < 1 1 -A 1 1 #2 [ I -B [ I jB 2 , so that tr(AB) is well-defined) and 

(A.6) det 2 , w (I w - AB) = det 2M {I H - BA) 

for all A G B(H,n'), B G B(H',H) such that BA G B 2 {W), AB G B 2 {%'). 
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